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Abstract 

We propose a new mechanism for magnetic field generation from inflation, by which strong 
magnetic fields can be generated on cosmological scales. These fields may be observable by cosmic 
microwave background radition measurements, and may have a dynamical impact on structure 
formation. The mechanism is based on the observation that a light nearly minimally coupled 
charged scalar may be responsible for the creation of a negative photon mass-squared (provided 
the scalar field coupling to the curvature scalar is small but negative), which in turn results in 
abundant photon production - and thus in growing magnetic fields - during inflation. 
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I. COSMOLOGICAL MAGNETIC FIELDS 

There is a growing evidence that magnetic fields permeate intergalactic medium [l|, 
0]. The galactic a — Q dynamo mechanism represents the standard explanation, according 
to which small seed fields in protogalaxies are magnified to a microgauss strength correlated 
on kiloparsec scales observed today in many galaxies. Yet the dynamo does poorly as 
regards to explaning the observed cluster fields, which are typically correlated over several 
kiloparsecs and can reach 10 microgauss strength. On larger scales the field strength drops 
quite dramatically, and it is characterized by a negative spectral index between —1.6 and 
—2 j^, consistent with the Kolmogorov spectrum of turbulence. 

The seed field that fuels the dynamo, which amplifies the galactic (and possibly cluster) 
fields, has either primordial origin, or it was created by the Biermann battery mechanism, 
operative at the time of structure and galaxy formation, when the Universe was about one 
billion years old B fl 3. The Biermann battery mechanism is operative in the presence 
of oblique shocks, which generate nonideal fiuid that violates proportionality between the 
gradients of pressure and energy density. Such a fiuid can produce the vorticity required for 
generation of seed fields, whose strength is typically of the order 10~^° gauss. While in the 
protogalactic environments these seed fields may be sufficiently strong to fuel the galactic 
dynamo mechanism, whether they can be used to explain the observed cluster fields is a 
highly controversial question. Moreover, these fields may be insufficientlystrong to explain 
the galactic fields observed in a couple of galaxies at higher redshifts |l|,y|- 

The quest for the fields correlated on much larger scales has so far been unsuccessful. 

n _ 

A notable exception is the result of j9|, where a local field of about 10 ® gauss correlated 
on megaparsec scale is quoted. If established, the existence of extragalactic fields would 
point at a primordial origin. Perhaps the most promising method for detecting primordial 
fields correlated over supragalactic scales is through their impact on the cosmic microwave 
background radiation (CMBR) jiol . 11 1, and possibly through their dynamical infiuence on 



large scale structure and galaxy formation. The current upper bounds from the CMBR are 
of the order 10~^ gauss, while the bounds from structure formation are somewhat weaker. 
Fields stronger than about 10~^ gauss at the time of galaxy formation (which correspond to 
about 10^^ gauss comoving field strength) are dynamically relevant. 

In this Letter we propose a novel mechanism for generation of large scale magnetic fields 
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during cosmological inflation. Our mechanism is based on a model with a rather minimal set 
of assumptions : We study the dynamics of (one-loop) quantum scalar electrodynamics 



oi assuraptions: 



in cosmological space-times, which can be naturally embedded into the 
Standard Model (where the role of the scalars is played by the charged components of the 
Higgs field), or into supersymmetric extensions of the Standard Model (the charged scalars 
are the sypersymmetric partners of the Standard Model fermions). 

II. SCENARIO 

We assume a standard model of the Universe, in which a period of inflation is followed by 
radiation and matter eras. The background space-time during inflation can be accurately 
approximated by the de Sitter metric, g^i, = a^rj^^ (r/^j, = diag[— 1, 1, 1, 1]) with the scale 
factor a = —l/{Hr]) (for t] < —r]H = ~l/-^)> where t] denotes conformal time, and H the 
inflationary Hubble parameter. After inflation the Universe undergoes a sudden transition 
to radiation domination, in which a = Hrj (for rj > rin). When inflation lasts a sufficiently 
long time, such that a de Sitter invariant solution is established, the dynamics of photons 
that couple to light scalar particles, which in turn couple nearly minimally to gravity, can 
be described by the Proca Lagrangian 

Cproc. = -\v'''v''''F>.uF,.- \aV^r^^'^A^A,+ 0(s°), (1) 
where = dfj,Atj — d^A^. The photon mass is given by 



where a = e'^/An is the fine structure constant, TZ = 12H^ is the curvature scalar of de 
Sitter space, m$ is the scalar mass and ^ describes the coupling of the scalar field to TZ. 
When < m| <^ and — 1/12 ^ ,^ < 0, both s and can be negative. We emphasise 
that there is nothing patological about this limit. Inde ed, the de Sitter invariant (Feynman) 
scalar propagator iAp[x,x') = Gpiv) va. D = A reads [12| 

where y{x;x') = aa' H"^ Ax'^ {x; x') denotes the de Sitter invariant length, and Ax'^{x;x') = 
~{\v ~ VI ~ '^^y + 11^ ~ ^ ~^ 0^) such that a negative small s introduces negative, but 
finite, correlations, which are at the heart of the mechanism discussed in this Letter. 



The Lagrangean ^ was derived in |l2'] by expanding in powers of |s| ^ 1 the one- 
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(order a) nonlocal effective action written in the Schwinger-Keldysh formalism. In 
proved that, when written in the generalized Lorentz gauge d^{a?ri^^ A^^ = at the order 
0{a, s^^), the effective action reduces to the Proca theory (|T|). 

Due to spatial homogeneity, the general solution for the photon field of the Proca equation 
can be conveniently written as a superposition of spatial plane waves A^{x) = ?7)e*'^'^, 
where = (^O)^- Then Eq is nondynamical and just traces the spatial components. The 
spatial componets correpond to the three physical degrees of freedom, and can be decom- 
posed into longitudinal e^^^ = k{k ■ e)/k'^ and transverse part ^ta- = ^Lfc- Fi'oi^ Eq. ((H), 
one finds the following equation of motion for the transverse polarization ^ during inflation 

{dl + fc2 + my ) e^ j:{r]) = , (4) 

where k = \\k\\. Demanding that the solution of this equation corresponds to a (circularly 
polarized) vacuum state for i] — oo, one finds Sj^ii]) = ^fc(^)(^g +^^|)/v^5 where 
(T = 1,2) are the transverse polarization vectors and 



A^, = li-nvM'\-kr^), z.= ^i-^, (5) 



Wronskian W[e.j.^,e*j^] = i. 



where Hi, is the Hankel function. The second solution is simply such that the 



III. VANISHING CONDUCTIVITY 

When conducitivity in radiation era is negligibly small, a smooth matching of the infla- 
tionary epoch solution (0) to the radiation era modes 

V2A; 



leads to the radiation era solution (c/. Ref. [nj). By matching Eq. Q to a linear combination 
of these solutions, one finds 

A = akAt\^=nH + l3lAz\^=r^^ , (7) 
A = afc(9,4)|,=,, + (3l{dr,Az)\^=^^ , (8) 



4 



where = Aj:\n=-r,„ and Aq = Solving for ak, (31 gives 

such that the gauge field evolves in radiation era as 

A! 

A^{r]) = Ao cos k{r] - r]H) + -j^ sin k{r] - r]H) {t] > ijh) ■ (10) 

For modes that are superhorizon at the end of inflation {k <^ H) one can use the small- 
argument expansion of the Hankel function Hl)^\k/H) to get (for z/ > 0) 

Ao = -i^^{2HY~h~'' + 0{k'',k"'+^) (11) 



A = {i^-1)hAo. (12) 

This result is used below to calculate the magnetic fleld spectrum in radiation era. Since 
matter era is a conformal space-time, the flelds (as well as their spectra) in matter era are 
inherited from the radiation era flelds ()10|). 



IV. LARGE CONDUCTIVITY 

An opposite extreme is a sudden increase of conductivity in radiation era, which occurs 
for example in the case of a rapid thermalisation after inflation, resulting in a large (thermal) 
conductivity. To leading logarithm (squared) in the coupling constant, the evolution of the 
photon fleld is then governed by the Bodeker-Langevin equation 

iaad, + e)e'^ir^) = a'eik,v), (13) 

where is deflned hj e^j: = J2t^'^^J- '^^^ stochastic force C,^ satisfles the following 
Markowian fluctuation-dissipation relation 

ieik, v)e'*ik', v'))~i2n)'5^niv - V')S^'^ (k - k') , (14) 

where T denotes the equivalent temperature (average energy of the plasma excitations). The 
factor a^, which multiplies .^'^ in Eq. (jl3|) . can be inferred from the fact that the transverse 
components of the current density are given by j = crE'^+C,'^, and from the transformation 
properties of the conductivity and the electric fleld. Eq. (fTSj) can be easily integrated, 

exp(^-J d7]"kyiaa))^dr]' + exp(^-J —dr]'je^ivH). (15) 



Using this and p4|) one can derive the equal time correlator 

+ exp(- f dr]'^^^^)e]:ivH)4'*iVH). (16) 

Assuming o", T oc 1/a (which is justified in and close to thermal equilibrium), one can perform 
the integrations, to obtain 



aT 



2k 
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{eUv)el'*m = -TT ^ - ^v) {2nfS^^'6^'\k-k') 



V aa 
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where Arj = rj — riH, such that the spectrum is neatly split into the thermal and primoridial 
contribution. 



V. MAGNETIC FIELD SPECTRUM 

When suitably averaged over a physical scale £ph = ai, the magnetic field operator can 
be defined as 

4(^,x) = J d^y\i{x~y,i)^{r],y), (18) 

where W denotes a coordinate space window function, which, e.g., can be chosen as W = 
[27r£2]-3/2exp(- p-y||V(2£2)), and 

Biv,y) = -,j ^,^k X X: e--aT-h.c.]. (19) 

The creation and annihilation operators dV^ and d'l obey the commutation relation, 

over the vacuum state, one arrives at 
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a^a^, 'J = (27r)^(5^^'5(^)(A; — k'), and a|'|0) = 0. When (fTH|l is squared and averaged 



{BKv,x)) ^ Ifm'Vs, 'PB = ^,^,\Adv)\', (20) 

where Vb = Vsijli ^) defines the magnetic field spectrum, and W(£, k) = f d^zV{z, rf) exp{ik- 
z) is the momentum space window function (W = exp(— /i:^£^/2) for the coordinate space 
window function W mentioned above). 



6 



When conductivity is vanishingly small, the spectrum in radiation era can be calculated 
by inserting Eqs. (fTUHT^ into ((201) • The result is (when a ^ 1) 



1 r\u)(2Hf'+\ 1,2 . 2rw 



(21) 



such that on superhorizon scales krj the spectrum (^1]) Vb oc k^~^'^, and on subhorizon 
scales the spectrum is oscillatory, with the envelope scaling as Vb oc k^'"^". This spectrum 
is shown in figure Q on a log-log plot. Note that scale invar iance is reached for z/flat = 5/2 
(superhorizon scales) z/flat = 3/2 (subhorizon scales). 
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FIG. 1: Magnetic field spectrum in radiation era (|21() for vanishing conductivity on a log-log scale 
for 1^ = 1 and z/ = 3/2. (Normalization is arbitrary.) 



In the limit of a large conductivity cr, the magnetic field spectrum can be split into the 
contribution from thermal excitations, and the primordial contribution. A new scale arises 
in Eq. p7|) . associated with the momentum 



kn 



(22) 



which we refer to as the conductivity scale. The thermal spectrum can be inferred from ()17p: 
In the ultraviolet cc k^ {k ^ k„), while in the infrared P]^ cc k^ {k <^ k„), implying that 
thermal excitations give rise to fields of negligible strength on cosmological scales (see jlT^). 
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The primordial contribution to the spectrum can be obtained from Eqs. ()17|) 
and PH). 



-expl — — a 



^4 27r3 ^'^-V Ha'^ V ^^^^ 

such that the spectrum is exponentially cut-off when k > k^. The flat spectrum is reached 
for i/flat = 5/2 for A; < k^^. This spectrum is shown in figure |2l on a log-log plot. Since typically 
in radiation era a ~ T ^ H/a (a ^ 1), fco-.ph = k^r/a ^ -f^(i) = H/a^, implying that a rapid 
growth in conductivity after inflation freezes out large scales magnetic fields and destroys the 
small scales fields (electric field spectrum is completely destroyed), such that no oscillations 
are present on subhorizon scales. Thus, the absence or presence of subhorizon oscillations, 
and observation of a conductivity scale, can be testing grounds for the conductivity history 
during radiation era. 




FIG. 2: Magnetic field spectrum in radiation l\2'6\i when a large conductivity sets in rapidly after 
inflation on a log-log scale for v = 1 and = 5/2. (Normalization is arbitrary.) 



VI. DISCUSSION 

The spectrum can be thought of as a function of the scalar coupling to gravity ^ and 
its mass m$, which are fundamental parameters characterizing a scalar field. To illustrate 
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this point more precisely, we show the spectral index n = 3 — 2z/ characterizing the spectrum 
envelope of Vb oc A;" in Eq. (j^D) on subhorizon scales, k/a > H{t) (or Be = {Bj)^/'^ oc 
(n ^ 0), Be oc -ln(£) (n = 0) in Eq. Between ^ = ^^rit = -ml/{12H^) and 

^ = ^3 = (a/Tr) — m|/(12if^), n = 3, equal to that of thermal spectrum. Above ^ > ^3, 
n drops, reaching an asymptotic value n rioo = 2 when ^ — 00. The spectral index on 
superhorizon scales is obtained simply by replacing n + 2^n = 5 — 2v. When compared 
with the vacuum spectrum, cc k"^ {k < H), a spectrum with n G (2,3) exhibits an 
enhancement of magnetic fields on subhorizon scales, which is due to a conversion of the 
electric energy (which is enhanced during inflation) into the magnetic energy during radiation 
era. The conversion is efficient provided radiation era is characterised by a low conductivity. 

0, Q 0, Q 



This mechanism was used in Refs. 
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Be oc can be obtained from inflation, which could be sufficient to seed the galactic 
dynamo mechanism. 

When C, G (~oo,^crit) however, gauge fields exhibit instability and are enhanced during 
inflation, such that the spectral index of subhorizon modes drops from n = 2 when ^ ^ — 00, 
to flat spectrum (n = 0), when 

Sir 12/72' 

to n < 0, when ^ G (^Hat, ^crit)- A negative spectral index n < —2 would imply a growth in 
magnetic field energy during inflation, resulting in a divergent magnetic (and electric) field 
energy, as it can be inferred from Eqs. (PT|). and with W = 1. A proper study 
of this case would require an inclusion of backreaction of the electromagnetic field on the 
background space-time, which is beyond the scope of this Letter. A spectral index n > —2 
on subhorizon scales (or > on superhorizon scales), results in an acceptable magnetic 
field spectrum on cosmological scales, whose dynamical impact on CMBR and large scale 
structure formation should be considered. 

Spectrum normalization today can be determined from magnetic field strength at a co- 
moving scale corresponding to the Hubble scale at the end of inflation {k = H). At the co- 
moving inflation scale ~ 3 m ~ 10~^^ parsec (corresponding to if ~ 10'^'' Hz ~ 10^'^ GeV) 
the field strength is about Bh ~ 10^^^ _ IQ^^^ gauss |20|. Therefore, if the spectrum is (al- 
most) flat, the magnetic field is potentially observable by the next generation of CMBR 
experiments (PLANCK satellite), as well as of importance for the dynamics of large scale 
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FIG. 3: Spectral index n for subhorizon modes in the low conductivity case (c/. Eq. (|21() ) as a 
function of the scalar coupling to gravity ^ (with a = 1/137 and the scalar mass m|, = 

structures of the Universe. 
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